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ABSTRACT 

We present a reliable estimate on the production rate of B c mesons in IS and 
25* states in the large transverse momentum region at hadronic colliders using heavy 
quark fragmentation functions derived within the framework of perturbative QCD. 
We also present the transverse momentum distribution for the B c mesons. The 
production rate is large enough for the B c mesons to be identified at the Tevatron. 
At the SSC/LHC the rate is so large that their properties can be studied in details. 

I. INTRODUCTION 

Heavy flavor production and decays are very useful for measuring CKM matrix elements 
of the standard model and for testing bound-state models for mesons and baryons. Ever 
since the first B meson was discovered, a lot of data have been coming out from, e.g., 
CLEO, ARGUS, LEP/SLC, FNAL, on the Bd and B u meson families. From the recent run 



at Tevatron the masses and other properties of the B s have been measured and confirmed 
The next family of B mesons will be the B c mesons made up of be. The 5 c -meson 
family differs from the J/ip and T families and from other B mesons because it is made 
up of a pair of heavy quark and antiquark of different flavors and masses. The J / ip and T 
families have played important roles in developing heavy quark bound-state models inspired 
by QCD. Being quarkonium systems of different flavors and masses be bound- states provide 
unique opportunities to test different bound-state models of QCD. The decays of B c mesons 
also provide rich sources to test the standard model, e.g. the measurement of |Vb c |, and 
enables us to see the interplay between strong and weak interactions. 

Since the physics of the B c mesons is so interesting, one would like to know how many can 
be produced in the present colliders (e.g. Tevatron) and in the future hadronic supercolliders 
(SSC/LHC). It is the purpose of this letter to present reliable estimate in the high p T 
region, by using the heavy quark fragmentation functions Di_^ Bc (z) which are based on 
perturbative QCD. We will summarize some features of these heavy quark fragmentations 
in Sec. II, and present the production rates and the pt distributions for the B c mesons in 
Sec. III. Finally we will conclude in Sec. IV. 

II. HEAVY QUARK FRAGMENTATION 

Previous estimates of B c meson production have based on perturbative QCD calculations 
for e + e~ colliders [[| and monte carlo studies for both e + e~ and hadronic colliders [|]]. The 
monte carlo estimates of the ratio <j{Bf) / a(bbX) are all of the order 10 -3 for LEP, Tevatron, 
SSC/LHC, and HERA. This fact leads us to think that the production mechanisms are all 
of the same nature. In the region of large p?, the major mechanism for producing B c or 
any other mesons is heavy quark fragmentation JJ,|| , in which a b antiquark is produced at 
large px by a hard-scattering process and it subsequently fragments into the meson. The 
differential cross section for direct production of the B c meson at high px can be factorized 
at leading order in ot s as 
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da(B c (p)) 



dz da(b(p/z, v))D- b ^ Bc (z, fi) , 



(1) 



where z is the longitudinal momentum fraction carried by the B c , and /i is a factorization 
scale. The physical interpretation is as follow: a heavy b antiquark is produced in a hard 
process with four-momentum p/z, and then it fragments into the B c meson with a longitudi- 
nal momentum fraction z. The fragmentation function D^ Bc (z) satisfies the Altarelli-Parisi 
evolution equation 



at leading order in a s . The factorization for B u , B^ and B s productions can be described 



fragmentation functions should be independent of the hard process by which the b is pro- 
duced. 

The fragmentation of b into B u , B d , and B s is a soft process, and can only be described 
by a phenomenological function ||. However, b — > B c requires the production of a cc pair 
and it is therefore a hard process which can be calculated using perturbative QCD @||. 
This perturbative QCD approach has been shown valid in calculating the fragmentation 
functions for heavy quarkonium productions ||, including the splitting of gluons and charm 
quarks into S- wave charmonium. The fragmentation functions Di^ Bc (z) derived in Ref. |2j 
need only the input parameters of a s , m^, m c , and the radial wave function R(0) of the 
bound state at the origin so that it has more predictive power. The initial fragmentation 
functions are given by 





in the same way as Eqs. ([I]) and (|2|), with the corresponding fragmentation functions. These 
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for the 1 Sq B c state, and 
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x [2 - 2(3 - 2r)z + 3(3 - 2r + Ar 2 )z 2 (4) 
-2(1 - r)(4- r + 2r 2 )z 3 + (1 - r) 2 (3 - 2r + 2r 2 )^ 4 ] , 

for the first excited state B* which has the spin-orbital quantum number 3 Si, and r = 
m c/ {fnb + m c)- Note that the scale in a s is set to be 2m c . The perturbative QCD calculation 
gives directly the fragmentation function at the scale /j,q = 2m c , which is the minimum 
virtuality of the gluon splitting into cc. Using the Heavy Quark Effective Theory methods 
@, it can be shown that the evolution up to the scale my, is trivial, which means that higher 
order radiative corrections will not introduce any logarithm of mb/m c . A convenient choice 
for the initial scale is Hq = mb + 2m c [[J, because the fragmentation functions for c — > B c , B* 
can then be obtained from Eqs. (|3|) and (|]) simply by interchanging mi, and m c . 

In the factorization scheme of Eq. ([]]), all the dependence on the momentum p is in the 
hard process a. Large logarithm of p/fi can be avoided by choosing the factorization scale 
to be of order p. The induced large logarithm of order fi/mi, in D(z) can be solved by 
evolving the Eq. (0). To leading order in a s only the b — > b contributes to the evolution. 
The fragmentation functions Di_> Bc (z) and Di^ B ,(z) at the initial scale [i and higher scale 
have been shown in Fig. 1 of Ref . H . 

III. RESULTS & DISCUSSIONS 

In this section, we use Eq. ([[]) to compute the direct production rates of B c and B* in 
hadronic collisions. Our calculation includes 

99 ~>bb, gb^ gb, and qq -> bb , (5) 

as the hard subprocesses for the inclusive production of the b. We choose the scale \i for 
the parton distribution functions and for a s to be the transverse mass of the 6, \jp 2 T i + rn 2 . 
We use the parameterization of HMRS (set b) [|8] for parton distribution functions. The 
running coupling constant a s (Q) is evaluated at 1-loop by evolving from the experimental 
value a s (rriz) = 0.12 ||, and given by 
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where bo = (33 — 2n/)/487r 2 , and n/ is the number of active flavors below the scale Q. 
The subprocess cross sections are convoluted with D(z,fi), as is in Eq. (jH). The functions 
D (z,fio) at the initial scale /xq are given in Eqs. (2) and (f|), and are evolved to the scale fi 
using Eq. (fj). For the initial fragmentation functions D(z,fi ) we use the input parameters 
of m b = 4.9 GeV, m c = 1.5 GeV, and \R(0)\ 2 = (1.18 GeV) 3 . 

The pt spectrum for the B c meson at Tevatron energy is shown in Fig. [I], with the 
acceptance cuts 

Pt(B c ) > 10 GeV and \y(B c )\ < 1 . (7) 

The corresponding spectrum for the B* is also shown on the same figure. The shape of the 
two spectra is very similar, because Di^ Bc {z) and Di^ B *(z) have similar shapes and differ 
primarily by an overall normalization difference of about 50%. The corresponding spectra 
at SSC(y / i = 40 TeV) and LHC(-\/i = 14 TeV) energies are shown in Fig. ||, but under 
slightly different acceptance requirements 

Pt(B c ) > 20 GeV and \y(B c )\ < 2.5 . (8) 

The integrated cross sections versus p^ m (B c ) are also shown in Fig. |^. The cross sections 
at the SSC are about three times as large as those at the LHC, and about two order of 
magnitudes larger than those at the Tevatron. 

So far we have only estimated the B c meson productions in 1 l So and 1 3 S\ states. Since 
the annihilation channel for the decay of the excited 5 c -meson states is suppressed relative to 
the electromagnetic and hadronic transitions, all the excited states (l 3 Si, 25*, IP, 2P, ID) 
below the BD threshold will decay to the ground state 1 1 5o by emitting photons or pions. 
They therefore all contribute to the inclusive production of B c mesons. A simple modification 
can be made to estimate the productions in 2 1 S and 2 3 Si states, by multiplying with 
the factor |i?2S , (0)/-Ris , (0) | 2 ~ 0.6 |10|]. Therefore, the curves in Figs. [I], |2], and |3] can be 



multiplied by 0.6 to get the productions for 2S states. To get the total inclusive production 
rate, however, we need to include P-wave and possibly D-wave contributions, which have 
not been calculated. Therefore, the cross sections presented here are rather conservative to 
the actual production rates. Table | shows the number of B c mesons that can be produced 
inclusively at Tevatron, SSC and LHC, including the contributions from IS and 2S states, 
with integrated luminosities of 0.025, 10, and 100 fb^ 1 respectively. It is also informative to 
give the ratio a(B c )/cr(bX), which is simply the fragmentation probability Jq 1 dzDi^ Bc (z). 
Adding the contributions from IS" and 2S states, the ratio a(B c )/a(bX) is about 1.5 x 10~ 3 , 
which is consistent with the monte carlo studies 

The factorization in Eq. ([!]) is correct up to the order (rub/pr) 2 , which explains why we 
impose a rather high cut on the B c mesons. In the low Pt region, mechanisms other than 
heavy quark fragmentation have to be taken into account. For example, production of pairs 
of bb and cc followed by recombination of b and c to form a B c meson contributes at low Pt 
region; whereas heavy quark fragmentation dominates at high p?- Similar conclusion can 
be found in the production of J/ip by heavy quark fragmentation [IT], which is dominant 
over the process gg — > ipg in the large pt region. There are other uncertainties arising from 
higher order QCD corrections, relativistic corrections of the bound state model, and the 
values of ot s and R(0) used. But the largest source of uncertainties comes from the values 
of rrib and m c employed because of the factor 1/mj? in the initial fragmentation functions as 
in Eqs. (|) and (§. 



IV. CONCLUSIONS 

We have used heavy quark fragmentation functions derived from perturbative QCD to 
calculate the production rates and the pt distributions for S*-wave B c mesons at Tevatron, 
SSC , and LHC energies. Imposing cuts of pt{B c ) > 10 GeV and \y(B c )\ < 1 at the Tevatron, 
and including the contributions from IS and 2S states, about 16000 B+ and 16000 B~ mesons 
should be produced for 25 pb _1 integrated luminosities. The corresponding numbers for SSC 
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and LHC with p T (B c ) > 20 GeV and \y(B c ) \ < 2.5 are 7.1 x 10 7 and 2.3 x 10 8 with integrated 
luminosities of 10 and 100 fbr 1 respectively. These B c mesons can be detected via the decays 
of the form J/ip + X, and in particular, via 

Bf -^i>^v t and Bf -> ^tt* (9) 

with if; — ■> . The first one has a distinct signature of three charged leptons coming off 
from the same secondary vertex and has a combined branching ratio of about 1%. The second 
decay channel has a smaller branching ratio of order 0.2-0.4%, but it has the advantage that 
the B c can be fully reconstructed. The production rate of B c mesons given above is large 
enough that it should be possible to use these decay modes to identify the B c mesons at the 
Tevatron, and to study its properties in details at the SSC and LHC. 
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TABLES 

TABLE I. Table showing the number of B c mesons that can be produced at Tevatron, SSC, 
and LHC, including the contributions from IS and 25 states, and under different p^ m (B c ) cuts. 
Another acceptance cut is |y(-B c )| < 1 (2.5) at Tevatron (SSC/LHC). The integrated luminosities 
are 25 pb _1 , 10 ftr 1 and 100 ftr 1 for the Tevatron, SSC, and LHC respectively. 
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FIGURES 

FIG. 1. The dependence of the differential cross section da/dpT(B c ) on the transverse momen- 
tum pt(S c ) for the ground state B c (l 1 So) and the first excited state B*(1 3 S±) at the Tevatron. 

FIG. 2. The dependence of the differential cross section da/dpT(B c ) on the transverse mo- 
mentum pt{B c ) for the ground state B c (1 1 Sq) and the first excited state B*(l 3 Si) at the SSC 
and LHC. 

FIG. 3. The integrated cross sections a (pt{Bc) > p^ m (Bc)) for the ground state B c (l 1 Sq) 
(solid) and the first excited state B*(l 3 Si) (dash) versus pf m (B c ) at the Tevatron, SSC, and LHC. 
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